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DIFFUSION OF CHARGED PARTICLES IN A RANDOM MAGNETIC FIELD 

JAMES A. EARL 

Department of Physics and Astronomy 
University of Maryland 
College Park, Maryland 20742 

ABSTRACT 

When charged particles move in a random magnetic field 
superimposed upon a relatively large constant field, their 
pitch angle distribution can be calculated to any desired pre- 
cision by an iterative approximation proceedure . Improved know- 
ledge of the pitch angle distribution and of the characteristic 
time for relaxation of anisotropy leads to an accurate expression 
for the coefficient of diffusion parallel to the mean field. 
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I . INTRODUCTION 

Scattering of charged particle trajectories by random magnetic fields 
controls the propogation of cosmic rays in the galaxy and in the solar system. 
The basic microscopic theory of this phenomenon is very complete (Jokipii, 1966; 
Hall and Sturrock, 1967; Hasselmann and Wibberenz, 1968), but the simplest 
application of this theory to diffusion has met difficulties that are not re- 
solved and the exploration of higher order transport phenomena has not yet 
begun. One manifestation of the difficulties is the fact that two different 
methods of calculating the coefficient of diffusion parallel to the average 
"dc" field give different answers. In the first method (Jokipii, 1966; 
Hasselmann and Wibberenz, 1970), diffusive anisotropies, treated as a small 
perturbation on the particle distribution function, are found by familiar 
proceedures similar to those used in quantum mechanics. The second method 
(Jokipii, 1968a), which involves an expansion of the distribution function as 
a series of Legendre polynomials, embodies a well known and often used pro- 
ceedure that is extensively discussed in standard references on transport theory 
Which method is correct? 

To answer this question, this paper develops transport theory in terms of 
eigenfunctions of the operator which describes pitch angle scattering (Section 
II) . Many aspects of this treatment which are closely analogous to the standard 
development in terms of Legendre polynomials will not be discussed in detail. 
Instead, specific references to the volume by Weinberg and Wigner (1958) will 
be given where appropriate . 

In section III, diffusion is discussed with special emphasis upon the 
important role played by the relaxation time for the lowest order anisotropic 
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eigenfunction. An exact evaluation of this parameter leads to a precise 
expression for the parallel diffusion coefficient. However, it will be seen that 
the diffusion coefficient obtained from perturbation theory represents a fairly 
satisfactory approximation accurate enough for many purposes. The derivation 
of the diffusion coefficient from a Legendre expansion leads to inaccuracies 
which put severe limits on the validity of this approach. Diffusion 

perpendicular to the average field is not affected by the considerations pre- 

t 

sented here* and, consequently, will not be discussed. 



II. EIGENFUNCTIONS OF THE SCATTERING OPERATOR 


Jokipii (1966) has demonstrated that pitch angle scattering can be 
described by a Fokker Planck coefficient of the form 


<p(u) 


- ^ - A IpI'-W) 


( 1 ) 


where p is the cosine of the pitch angle. The parameter A is given by 

' ' Vk q 

4 ■ 75 <-f-> < 2 > 

R o 

where to is the relativistic Larmor frequency, and V and R are, respectively, 
particle velocity and rigidity. The energy density contained in fluctuations 
of magnetic-field components perpendicular to the average field of wave number 
k in the interval dk is assumed to be Q (k /k)^ dk so that Q represents the 

XX O XX 

spectral density at the reference wavenumber k Q . In an interplanetary context 
where spatial inhomogeneities convected by the solar wind past a satellite 
observing platform give rise to temporal fluctuations of magnetic intensity 
characterized by a power spectrum P(f) = P^Cfo/f) , the parameter A can be 


expressed as, 


V V 2rtVf q 

A = -2- P ( °) 

_2 xxV 0) ' 

R wo 


(3) 


where is the velocity of the solar wind, and P^ is the spectral density 
of perpendicular magnetic field components at the reference frequency f 
(Jokipii, 1967; Jokipii and Coleman, 1968, Jokipii, 1971). 

The Boltzmann equation can be expressed in terms of the Fokker- 


Planck coefficient cp(p) as 
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where f(p,z,t) is the distribution function, z is distance parallel to the mean 
field, and t is time. In the treatment that follows, the streaming term 
pV3f/dz will be regarded as a perturbation which is usually but not necessarily 
small . The equation 


5 5R K 2 

5? ” ST + <tt> *K ■ 0 

K. 


(5) 


defines eigenfunctions R^p) and eigenvalues (2 /t^) which are useful in describing 

the time evolution of the distribution function. Physically, T is the 

K 

relaxation time required for an anisotropy proportional to to decay by a 
factor of (1/e). Sturm Liouville theory ensures that the R^ exist and form 
an orthogonal set satisfying boundary conditions that the finite value 
R^(l) = ± R^(-l) and that R^CO) is finite. 

In the spirit outlined above, we write Equation 4 as 


if 1 5 m 3£. , , 

2 55” 


( 6 ) 


where e is the perturbation and express both e and f as eigenfunction series 


f = Ef K (t,z)R K (p) 


(7) 


where the expression 


e = Ee K (t,z)R K (p) 


( 8 ) 


+1 

f. e 


\ d)1 


e„ = 


K +1 


(9) 


/ V 

-1 


dp 


which follows from the orthogonality of the R^ describes the coefficient e^. 

A similar expression gives the coefficients f^. Equations (7) and (8) can be 
inserted in equation (6) to yield, upon invoking the orthogonality of the R^,, 



6 


a set of independent equations for the coefficients f 


K 


Sf K 1 

at + ( t k ) f K e K 


( 10 ) 


whose solution is 


f K (t,z) = f K (t=0)e 


./U 


z)e 


-(t-s)/T 


K 


ds 


( 11 ) 


Equation (11) shows that the value of f at time t is primarily dependent upon 

K 

values of e within a period of width T preceeding t. When t » T and when 

K. K 


does not change significantly in a time T v , (l/e)(Se/St) « (1 /t^.), then f. 


reaches an equilibrium value 


K 


K 


K 


f = e T 
K K K 


( 12 ) 


Under these conditions, equation (12), together with general theorems which 

2 

require that the eigenvalues (2/t ) increase with K at least rapidly as K , 

K 

guarantee that the series for f (Equation 7) converges absolutely and uni- 
formly (Courant and Hilbert, 1953, p. 427).' By the same token, the approach 
to equilibrium of higher order eigenfunctions is more rapid than of lower 
order ones. (See Equation 11.) 

In the case of isotropic scattering, f = A(l-y 2 ), the eigenfunctions 
are Legendre polynomials P (y) with eigenvalues (2/r ) = K(K+1)A. When f is 

K lx 

complicated, the are found by standard methods which also yield useful esti- 
mates of the eigenvalues. The qualitative behavior of the is defined by 
theorems which also apply to the P . Consequently, many features of 
the can be predicted by analogy with the P^. For example, eigenfunctions 
having even indices are even functions of y while those having odd indices are 
odd functions, the number of zeros between y = -1 and y = +1 is equal to the 
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index, etc. The lowest eigenfunction R = const, describes the isotropic com- 

o 

ponent and is, of course, identical to P . The fact that the relaxation time 

o 

is infinite means that an isotropic distribution is stable. 

Because the ratio of successive terms in equation (7) depends not 
only upon the index K but also upon the boundary conditions and coefficients 
e^, no general statements can be made about how rapidly the series for f converges. 
However, for the case considered here, e a y, the properties of R^ just described 


imply that, for K > 1, the e decrease rapidly with increasing K but are not, 

K 


necessarily , identically zero as they are when R^ = P . This fact 


ensures 


the 


validity of the basic assumption of the diffusion approximation^that terms 
above R^ can be neglected. In contrast, when the diffusion problem is described 
by another arbitrary set of orthogonal functions, f is essentially represented 
by the isotropic component plus a series expansion of R^ in terms of the new 
basis set. In general, the convergence of this series, which depends only 
upon the nature of R-^ and not upon physical parameters, is not as rapid as 
that of the series in R^. An example that will be discussed in detail later 
is the Legendre series for R. when q ^ 2 in which the coefficients of P 
decrease as (1/K) , much less rapidly than the (1/K 2 ) dependence implied by 
equation (12). These facts single out the R^ as a preferred set of orthogonal 
functions whose use leads to greatly simplied analysis. 



8 


To approximate Rp we invoke the iterative method of Stodolla and 
Vianello (Hildebrand, 1949, Chapter 5). In the first iteration, a trial 
function w satisfying the boundary conditions but otherwise chosen arbitrarily 
is inserted in Equation 5 


3 dx 2 . v 

55 <p 55 ' ' ? " (p) 


(13) 


to obtain, after two integrations, an improved estimate x of the eigenfunction 


x 


h* * 

(M) = ( T } / f W(P) dP 


(14) 


0 v 

In the first integration over p, the constant must be chosen so that, at 
v = +1, cp(v) (dx/3v) = 0. (See Equation 1.) To satisfy a similar restriction 
at v = -1, w must be an odd function of p. In the second integration over 


V, the constant is chosen to make x an odd function of p. An accurate 
(Hildebrand, 1949, Chapter 5, Eq. 89a) 
estimate/ for the eigenvalue (2/t) is 


+1 


(2 )m zL 
V +i 


y w 2 dp 


(15) 


$ gw dp 

-1 

where the function g - (t/2)x is the double integral which appears in equation 
(14). This process can be repeated by using x as the trial function, replacing 
w, to calculate an improved estimate y of the eigenfunction. 

If w = p and cp is described by equation (1) with q < 2, equation (14) 
gives: „ ..2 % 


.(u) = (2) I r dvLl - 
V 2 J cp( v) 

2-9 


0 


= &) 


TA' 2(2-q) 


(16) 
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an expression identical to that obtained by Hasselman and Wibberenz (1970) 
from perturbation theory. In equation (16), as well as in all other expressions 
for the first eigenfunction presented in this paper, only positive values of 
p are considered. Since we are dealing with odd functions, corresponding 
results for negative p are easily obtained by reflection. The approximate 
eigenvalue calculated from equation (15) is 

(fj[> = | <2-q)(4-q) (17) 

from which the diffusion coefficient first derived by Jokipii (1966) can be 
obtained by the procedures outlined in the next section. 

In the second iteration, using equation (16) as the trial function, 
equation (13) is satisfied identically for p -* 0 provided that 

( ?^A ) = C 2 -q)( 3 -q) (18) 


which is, therefore, an exact expression for the eigenvalue, 
approximate eigenfunction, accurate enough for most purposes, 

/ 1-q 1 - V 3 ~ q 


y (p) = c i J v 


dv 


l - 


The corresponding 
is given by 

(19) 


where the normalization constant 

(1/C i ) = 2 [[♦( 5 r a -) - ♦c 1 ?) J (2°) 

can be defined in terms of the digamma function t|i in such a way 

that y(l) = 1. In the limit p -* 0, 

y(y) = T2^T * 2-q (21) 

while an accurate expression when p -* 1 is 

3-w - 1 -■ r ^ ] < 22 > 
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In Figure 1, the function y is plotted vs. p over a range from q = 2, 
where the eigenfunction reduces to a step at the origin from -1 to +1, through 
q = 1, where the Legendre polynomial = yt is obtained, to q = -1, where a 
pronounced peaking in the forward direction (p=l) is evident. For q > 2, the 
treatment given here breaks down and the eigenvalues diverge but for all other 
values of q including negative ones the method converges to well behaved 
eigenfunctions. For q < 1, the eigenfunction has zero slope at the origin but 
this slope is infinite when 1 < q < 2. In the latter range, the shape of the 
eigenfunctionis evidently quite sensitive to the parameter q. 

It is appropriate in the present context to introduce the Legendre 
expansion of the eigenfunctions: 


V") * ] Vj 00 <23) 

and to express y(p) in terms of its coefficients a^ . This representation not 
only provides a useful method of evaluating certain parameters discussed later 
but also it makes explicit the relationship between the approach given here 
and the classical approach based upon Legendre expansions of the distribution 


function. When equation (19) is inserted into the 


formula for 


Legendre coefficients and when the order of integration over p and v is 


interchanged, the coefficients are given by, 

1 . .. ..3-q 1 

r 1-q 1-V H r 


a lj = C 1 (2j+1) / dV V 


i-v ^ r 

i - 2 i 


Pj(n) dp 


' 1 j(J+l) 


L " q < i - v3 ' q > fcy v > 


where the second integral, which involves only powers of V, leads to the following 
expressions for the first few coefficients: 
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3 < 3 -q> c i 

a ll 2 (2-q) (5-2q) 

7 (i-q)(3-q) (ll-4q) ^ 
a 13 = 4 (2-q) <4-q) (5-2q)(7-2q) 

n (q-l)(3-q)(8-3q)(71+138q-36q 2 )C 1 
a 15 = S4 (2-q) (4-q) (6-q) (5 — 2q ) (7-2q) (9-2q) (25) 

The coefficients for even indices are zero. The expressions in equation (22) 

are such that no singularities occur for q < 2.5. (Note that, as q *• 2, C 1 (2-q)-) 

When q = 1, a^ = a^,. = 0, and the eigenfunction reduces to P^(p). When 

q = -.4588, a^ changes sign. 

In Table I, numerical values of C^, given by equation (20), of 
C^/ (2-q) and of a^, a^ and a^,., given by equation (25), are 

tabulated for several values of q. The fact that C^/(2-q) deviates from one by 
no more than 10$ demonstrates that equation (21) is a surprisingly accurate 
representation of the eigenfunction. An equivalent statement is that even 
the first approximation described by equation (16) is fairly accurate. In 
Figure 2, which presents values of y given by equations (21) (dotted line) and 
(22) (dashed line) for q = 1.9, the Legendre series including terms up to 
(solid line) displays large oscillations around the correct eigenfunction. This 
manifestation of the Gibbs phenomenon is to be expected in a range of q where 
the eigenfunctions are virtually equivalent to the step function reached in the 
limit q -* 2. In fact the coefficients for q = 1.9 are nearly equal to those of 
the Legendre expansion of a step function: a^ = (3/2), a^ = -(7/8) and 
& - (11/16). In contrast, at q = 1,5 (Figure 3), the Legendre representation 
is relatively accurate. The accuracy of the Legendre expansion has an important 
bearing on the applicability of classical transport theory, a point that will 
be discussed further in the next section. 
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III. DIFFUSION EQUATIONS 

When is finite, the approximate solution of the Boltzmann equation 
obtained by keeping only the R Q and terms in equations (7) and (8) and by 
neglecting the time derivative (df^/dt) is the diffusion equation. (See 
Weinberg and Wigner, 1958, p.231). In this case, equation fl.0) gives: 


3f 3f 

2. — e — -u i. 

3t 0 " V 01 3z 


3f 

a! 




- ' f l - e i = “ V 01 rf* 


(26) 


(27) 


where the quantity Vq^ is a "matrix element" in the eigenfunction expansion of 
the perturbing term, 

„ a f V -a tt K 

e = -Vn ^ Rj V KJ (28) 

r+i 


I- 




rr +i 2 

f+i 

V- 

V 

J-l 

J-l 


n4 


V = V 
KJ JK 


Note that, because p is an odd function, 


(29) 


= 0, while 


V K,K+F ^K,K+3 etc are finite except for unusual cases when pR^ happens to be 
orthogonal to R^ + ^, etc. This coupling, through p, between odd and even eigen- 
functions is a ubiquitous feature of transport theory (Weinberg and Wigner, 1968, 
p.246). The diffusion equation results from inserting f^ calculated from 
equation (27), neglecting (3f^/3t), into equation (26) 


3f 


3f 


_ a 2 Sf o a . 

3t ~ ^ 1 V 01 3z ~ 3z 3z 


(30) 


with the diffusion coefficient D given by 
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K 


D - Voi ■ *V* L St 




r 1 2 

L v* 


(31) 


where the fact that R^ = 1 has been used in applying equation (29) . The inte- 
grals in equation (31) can be expressed in terms of the Legendre coefficients 
given by equation (25) 


f 


+1 


-1 


uRjdp = (2/3 a n 


(1 /C 2 )2 = ^ R 2 dp = (2/3) a 2 n + (2/7) a^ 3 + (2/11) aj 5 +. 


(32) 


where is evidently the constant multiplier required to convert R^ into a 
normalized eigenfunction. Numerical values of C 2 calculated from equation (32) 
appear in Table I. The diffusion coefficient is given by; 

D=(1/j^ 1 V 2 [3(V Q1 /V) 2 3 rj (1/3)>' 1 V 2 (33) 

where the second approximate equality follows from the fact that the quantity 

3(V 0 i/V) 2 = [f(a u ) 2 (C 2 ) 2 Ml +f (a l3 /a ll )2 + ll (a 15 /a ll )2 + " (34) 

is virtually equal to one. (See Table II.) 

Equations (33) and (34) make evident the fact that D is critically de- 
pendent upon but is quite insensitive to the exact shape of the eigenfunction. 
When t^v is defined as the mean free path \ the elementary formula D ~W/3 is 
recovered from equation (33). 

With the exact value of t given by equation (18) the diffusion coefficient 
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calculated from equation (33) is, 

y 2 2 £ 3 (V Q1 /V) 1 _ (2/3) vf 

" A 3 (2-q) (3-q) ~ (2-q)(3-q) A 


(35) 


an expression which can be compared to the result calculated with given by 
equation (17) (This result is identical to that obtained from the perturbation 
approach discussed by Jokipii (1966) and by Hasselman and Wibberenz (1970).) 


D 

P 


1 

(2-q) (4-q) 



(36) 


and to the expression derived from a Legendre expansion neglecting all terms 
beyond (Jokipii, 1968) (Weinberg and Wigner, 1958, p.233), 


D 


L 


= q(q± 2 ). v! 

9 A 


(37) 


To illustrate the relationship between these three expressions for the 

diffusion coefficient. Table II gives the ratios (D / D) and(D /D) as a function 

P L 

2 

of q. Also given in Table II is the quantity [ 3(V^/V) ] which is the ratio 

between the exact and approximate forms of equation (35) . The exact value of D 

+ for 1 < q < 2 

is a little smaller than that given by equation (36) /but deviates by no more 
than 25$ for 0<q<2. This good accuracy reflects the fact that equation 

(36) gives the first approximation in a rapidly converging sequence. The 
result derived from a Legendre expansion is not very useful since equation (37) 
does not predict the divergence at q = 2 and incorrectly gives zero diffusion 
for q = 0. These discrepancies occur because the Legendre derivation arbi- 
trarily sets equal to zero the coefficients etc • When these terms 

are taken into consideration, the diffusion coefficient can be represented by 
an infinite series in which equation (37) appears as the first term but this 



term is not necessarily the largest one. On the other hand, when q = 1, 
equations (35), (36) and (37) give identical values of D. In this special 
case, all three methods of calculating the distribution function are exact. 
The streaming flux S is given by the integral 

C +1 

S = (1/2)V ^ pf dp (38) 


which reduces to 


S = (l/3)a 11 Vf 1 + (1/3) a 3] Vf 3 + ... 


(39) 


an expression which involves only the coefficients a^ of the first Legendre 
polynomial = p. In contrast to the situation in classical transport theory 
where p is orthogonal to all other Legendre functions, all odd contribute 
to the flux. However, for reasons given earlier, the terms above f^ will 
almost always be negligible. 

The use of equations (27), (29) and (31) shows that the familiar expres- 
sion for the diffusive flux. 


= -n 


S = -D 


(40) 


retains its rigorous validity. On the other hand, an often used relationship 
between flux and anisotropy 6 which is based upon a Legendre expansion neglecting 


terns above 


6 

L 


max min 

f +f . 
max mm 



( 41 ) 


must be replaced by the correct expression 
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valid when f^»f , f,_, etc. For q = 1 equations (41) and (42) are identical 
but they give different results for q 4 1. For q ~ 2, the anisotropy given 
by equation (42) for a given flux is 5C$ smaller than would be expected from 
equation (41). (See the numerical values of a^ appearing in Table I.) These 
remarks apply to diffusive anisotropies only. Convective anisotropies, such as 
those discussed by Forman (1970), result from a term in the distribution function 
accurately proportional to p and, consequently, are described by equation (41). 
These difficulties could be avoided by defining the anisotropy in terms of the 
coefficient of . This definition lacks the direct relationship with observed 
intensities embodied in the first equality of equation (41) but it is entirely 
workable and, in fact, corresponds to what is measured by the usual Fourier 
analysis of periodic variations in the rate recorded by a rotating detector. 

When equations (26) and (27) are solved retaining the derivative (df^/dt) 
and assuming that does not depend upon z, we obtain 

a^f . a z f , af 

^ (V 0 i) at z d at ^ ' 

a version of the well known "Telegraphers equation" which embodies both diffu- 
sive phenomena characterized by D and wavelike propogation at velocity Vq^ . 
Equation (43) is analogous to a result of classical transport theory (Weinberg 
and Wigner, 1958, p.235) which has been discussed in the context of inter- 
planetary propogation of solar cosmic rays by Axford (1965) and by Fibich and 
Abraham (1965) . These treatments obtain for the propogation velocity (VA/3) 
which coincides with Vq^ only for q = 1. However, the deviations of Vq^ from 
(V/V3) are quite small as can be seen from the values of (a/^Vq^/V) given in 
Table II. Note that, for q ^ 1, the velocity is always a little smaller 
than (V/ . 
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IV DISCUSSION 

The above analysis has treated in detail only the special case of a 
power law spectrum of magnetic irregularities described by the two parameters 
A and q. More complicated cases involving additional parameters could also 
be analyzed by the methods presented here. However, it is evident from 
equations (14) and (15) that both the first eigenfunction and the relaxation 
time are given by integrals over p which are rather insensitive to the exact 
behavior of cp. This insensitivity implies that the above results can be used 
with confidence in describing the interplanetary propagation of 

cosmic rays whose diffusion is controlled by a portion of the 

irregularity spectrum characterized by a power law with index q ~ 1.5 (Jokipii 
and Coleman, 1968). On the other hand, caution should be exercised in applying 
these results to situations in which propagation is controlled by a 

steeper power law (Sari and Ness, 1969) with an index close to the 

critical value q = 2 where divergences appear. Hasselman and Wibberenz (1970) 
have suggested that, for q > 2, the diffusion picture breaks down and convective 
effects dominate. The observed steepening of the interplanetary magnetic power 
spectrum implies that diffusive anisotropies at high energies are more 
concentrated around the average field direction than at low energies. Because 
the eigenfunction is quite sensitive to the value of q (See Figure 1.), this 
difference might be observable. 

Jokipii (1968b) has already pointed out that the small amount of 
scattering near p = 0 implied by a Fokker-Planck coefficient described by 
equation (1) leads to a much larger relaxation time than that predicted by 
setting equal to the collision time as is suggested by a result of elementary 
kinetic theory. Equation (18), which predicts that -* 00 as q -* 2, provides 
a rigorous basis for these ideas . 



18 


The derivation of a nearly exact expression (Equation (35).) removes 

any conceptual uncertainty about the value of the diffusion coefficient but 

the numerical results are not very different from the first order approximation. 

(Equation (36).). The following formula, based upon equations (3) and (35), 

2 

gives the parallel diffusion coefficient D in cm /sec: 

2 -, 


D = 


2.22x10 


28 


(2-q) (3-q) 


3 ]( - B — T 

L -1\2. 09x10 


q v q-i 

W 


XX 


PR 


2-q 


(44) 


where p = (V/c), R is particle rigidity in GV, V is solar wind velocity in 

w 

km/sec, B is the average field in gamma, and is the spectral density in 
2 

(gamma) /Hz of perpendicular fluctuations at a reference frequency of f q Hz . 

Interpolation between the values given in Table II will yield an adequate 

2 

estimate of 3 (Vq^/V) . This expression is comparable to that given by Jokipii 
(1971, Equation (62).). In using equation (44), note that P^ x refers to a 
frequency spectrum defined, in terms of the Fourier transform of an auto- 
correlation function, for negative as well as positive frequencies. Observations 
are often given with the contributions of positive and negative frequencies 
combined. In this case, the reported spectrum is 2 P . 

Of greater significance, perhaps, than any of the above points is that 
the discussion of eigenfunctions given here provides a starting point for the 
investigation of higher order diffusive effects and of phenomena encountered 
when q > 2. 

This paper was written while the author was on leave at the California 
Institute of Technology where the hospitality and criticism of J. R. Jokipii 
were much appreciated. This work was supported, in part, by the National 
Aeronautics and Space Administration under Grant NGR 21-002-066. 
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FIGURE CAPTIONS 


Figure 1; The shape of the function y, obtained by numerical integration of 
equation (19), is strongly dependent upon the parameter q. 

Figure 2; For q = 1.9, the function y is accurately given near pi = 1 by equation 
(22) (dashed line) and near >u = 0 by equation (21) (dotted line) . The Legendre 
expansion (solid line), including terms up to P^, is not a very accurate re- 
presentation of y. 

Figure 3: For q = 1.5, all three representations of y give fairly accurate 


results . 
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Table I 

Numerical values of some parameters 


q 

C 1 

C 2 

C^y (2-q) 

a ll 

a l3 

a l5 

0 

1.562 

1.668 

.781 

.703 

.322 

-.036 

0.5 

1.326 

1.458 

.885 

.829 

.207 

-.011 

1.0 

1.000 

1.225 

1.000 

1.000 

.000 

.000 

1.5 

.548 

.972 

1.096 

1 .233 

-.359 

.213 

1.8 

.217 

.826 

1.086 

1.397 

-.662 

.471 

1.9 

.105 

.782 

1.055 

1.450 

-.771 

.577 

2.0 


.744 

1.000 

1.500 

-.875 

.687 

Equation 

(20) 

(32) 

(21) 

(25) 

(25) 

(25) 



Table II 

Comparison of Diffusion Parameters 


q 

( D p / D ) 

<V D > 

[ 3 ( V 01 / V ) 

0 

1.23 

0 

.917 

0.5 

1.10 

.80 

.974 

1.0 

1.00 

1.00 

1.000 

1.5 

.94 

.68 

.957 

1.8 

.92 

.31 

.887 

1.9 

.91 

.16 

.857 

2.0 

.90 

0 

.831 
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Figure 2 





Figure 3 



